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Abstract 

Among the several proofs known for X^^il/"^ = ""^/S, the one by Beukers, 
Calabi, and Kolk involves the evaluation of j^j^ — x^y^)dxdy. It starts 
by showing that this double integral is equivalent to | X^^iV"-^; ^-^^d then a 
non-trivial trigonometric change of variables is applied which transforms that 
integral into JJ^, 1 dudv, where T is a triangular domain whose area is sim- 
ply 7r^/8. Here in this note, I introduce a hyperbolic version of this change 
of variables and, by applying it to the above integral, I find exact closed-form 
expressions for [sinh~^ (coshu) — u] du, [u — cosh~^ (sinhit)] du, and 
/c^2 In(tanhu) du, where a = sinh^^(l). From the latter integral, I also derive 
a two-term dilogarithm identity. 
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1. Introduction 

The Riemann zeta function ^(s) is defined, for complex values of s, 5ft(s) > 1, 
by C{s) ■= J2^=i ^/n^- For integer values of s, s > 1, the first zeta value is 
C(2) — X^j^Li l/"-^- Euler was the first to estimate the numerical value of this 
series to more than 5 decimal places (1735), as well as to determine an exact 
closed-form expression for it (1740), as given by jl[: 

Lemma 1 (Euler's result for C(2)). = 7^2/6. 

From Euler's time onward, this result has been proved in several forms, from 
elementary to complex ones pj, Among these proofs, at least two involve the 
evaluation of a double integral over the unit square and these are the ones in 
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which we are interested here. Let us briefly describe these two proofs for a 
better exposition of the hyperbolic approach and its consequences. 

Proof (Apostol's proof). The short proof by Apostol (1983) can be found 
in Ref. (j). He starts by expanding the integrand of Bcukers' integral / := 
/oX^ i-xy dxdy as a geometric series, in order to show that it is equivalent to 

I 2 2 

^ n—l 

the interchange of integrals and sums being justifiable by the fact that each inte- 
grand is nonnegative and the sums converge absolutely (see the only corollary of 
Theorem 7.16 in Ref. By applying a simple linear change of variables cor- 

responding to the rotation of the coordinate axes through the angle 7r/4 radians, 
as given bjH 

U — V U + V 

X = — ■=- and y = — ■=- , 

directly to the Beukers' integral, one finds [1, @: 

pV2/2 arctanf , ) /-v^ arctan(^^=%) 
1 = 4 + 4 A^lldu. 

Jo \/2-v? J v2--u^ 

Now, the substitution u = \/2 sin6', followed by the use of the trigonometric 
identity arctan(sec6' — tan6') = 7r/4 — 6/2 on the second integral, yields / = 
7rVl8+ 7rV9. Therefore, Y.n=i V"^ = 7rVl8+ = tt'^/Q. □ 

The other proof, in which we are more interested here, was given by Beukers, 
Calabi, and Kolk (1993) 

Proof (BCK proof) . Similarly to Apostol's proof, Beukers, Calabi, and Kolk 
(BCK) start by showing that 

^ ■= 1 ^dxdy= / / Y.^'"y^"dxdy 

JoJo ^ ^ y Jo Jo 

00 „1 ^1 00 ^ 

Again, as each integrand is nonnegative the absolute convergence allows the 
interchange of the sum and the integral. Since J2^=i V(2"-)^ — j J2'^=i V"-^: 



^Note that 1/(1 — xy) is a positive real number at every point in the unit square, the only 
exception being the point (1, 1), where it is undefined. Of course, usual methods for treating 
improper integrals apply to this case. 

^For this change of variables, it is easy to see that \ J\ = 1, J being the Jacobi determinant 
d(x,y) _ dx/du dx/dv 

~ dy/du dy/dv ' 
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then X^S^o (2n+i)^ ^ I TT?- They then evaluate the unit square integral 

that defines K by applying the following non-trivial, trigonometric change of 
variables Q: 

sin u sin v 
X = , y = . 

cos V COS u 

For this change of variables, one easily finds that |J| = 1 — tan^ u tan^ w = 
1 — x^y^, which reduces the integral for K to jjj, 1 du dv where T is the domain 
in the uu-plane defined by {(u,v) : u>0, v>0, u + v< tt/2} 0, 0, |3|- 
This double integral exactly evaluates to 7r^/8, i.e. the area of the triangle 
corresponding to T, which implies Lemma[T] □ 

Here in this work, I introduce a hyperbolic version of BCK change of vari- 
ables for double integrals over the unit square. By applying this new change 
of variables to /q/q 1/(1 — x'^y'^)dxdy, I derive exact closed-form expressions 
for some definite integrals. From one of these integrals, I deduce a closed-form 
expression for Li2(V2-l)+Li2(l- V2/2). 



2. A hyperbolic version of BCK change of variables 



Let us introduce the following hyperbolic version of BCK change of variables 
for double integrals over the unit square: 



sinh u sinh v 

X = ; — and y 



coshw 

The corresponding Jacobi determinant is: 

coshu sinh It sinhw 



coshi 



(2) 



J 



coshw 

sinhu sinhw 
cosh^ u 



cosh V 
coshu 



coshu 



1 — tanh^ u tanh^ v. 



As < a;?/ < 1 in the unit square, then 1 — x^ y"^ > 0, then | J| is also equal to 



1 — tanh^ u tanh^ v. Therefore 



1 /.I 



y) dx dy — 



F{x{u,v),y{u,v)) \J{u,v)\ dudv 

F{u,v) (l — tanh^ u tanh^ u) dudv, (3) 



where S is the cusped hyperbolic 'quadrilateral' on the right-hand side of Fig.[TJ 
We are assuming that this change of variables is a diffeomorphism, which is 
reasonable since the functions x = x{u,v) and y = y{u,v) given in Eq. as 
well as 



u = u{x, y) = cosh 



1 — a;^ y2 
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and 

V = v{x, y) = cosh" 



1 i + y^ 



2-2 y2 



which define the inverse change^ are continuously difi'erentiable in all points 
of their domains (see the hachured regions in Fig. [T]) , the only exception being 
point {x,y) = (1, 1). As indicated there in Fig. [H S is bounded by the axes u 
and V, the curve v — f{u) = sinh~^ (coshu), above the diagonal v = and the 
curve V — g{u) — cosh^^ (sinhu), below the diagonal. The parameter a can be 
readily evaluated by noting that it is the abscissa of the point where v ~ g{u) 
intersects the u-axis. There in that point, one has cosh~"'^(sinha) — 0, thus 
sinha = 1. Since sinh""'^ x = In (x + Va;^ + 1 ), Vx G M, then 

a = sinh"^ (1) = ln(l + Vi). (4) 



3. New definite integrals and a dilogarithm identity 

Let us now apply our hyperbolic change of variables to the unit square 
integral as defined in Eq. ([1]). 

Theorem 1 (First definite integral). 

/ Tsinh""'^ (coshu) — ul du = — . 
■/o ^16 

Proof. From the BCK proof, one knows that ^^^^ 1/ {1 ~ x^y"^) dx dy ~ 7r^/8. 
By combining this result with (|3]), one has 

75 T^ — f 1 — tanh^ u tanh^ v] du dv 



si — tanh u tanh v 
Idudv = area of S. 

S 



The area of the domain S can be equalized to a single definite integral by noting 
that S is symmetric with respect to the diagonal v = u (see Fig. [1}. Then, by 
taking only the upper half of S into account, one has 

area of 5 = 2/ / Idvdu = 2 i [sinh~^ (coshu) — u] du , 
Jo Ju Jo 

which leads to the desired result. □ 



•^Note that the real function sinh^-'^a; is the inverse of the bijective function sinhx for 
all a; £ R, as usual. However, since the function cosh a; is not bijective over the entire real 
domain, we can reduce its domain to the nonnegative reals and then define its inverse as 
cosh" ^ X : = In ( a; + x'^ — 1 ) , V x > 1 . 
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Theorem 2 (Second definite integral). 

J [u~ cosh"^ (sinhu)] du ^ — - - In^ (^1 + V2j . 

Proof. Similarly to the previous proof, by taking the lower half of S into 
account and noting that it is bounded by the w-axis, the diagonal v = u, and 
the curve v — g{u), as seen in Fig.[Tl one has, from Eq. that: 



area of S' = 2 



a 
Y 



1 dv du 



la J g(u) 

Again, since the area of S is 7r^/8, one finds 



f oo 

/ [u — cosh^^ (sinh u)] du . 

J a 



\u — cosh ^ (sinhw)] du 



1 /tt 



2 V 8 



16 



□ 



Theorem 3 (Third definite integral). 

[ In(tanhz) = i In^fl + \/2') - - 

Ja/2 ' ' 4 V J h 



2 

16 



Proof. From Theoremfl] we know that the area of S evaluates to 



-'it 



/(«) 



I dv du — — 



where f{u) — sinh~^ (coshu). Let us rotate the coordinate axes by the angle 
7r/4 rad, as in Apostol's proof. Since "pure" rotations about the origin always 
yield an unitary Jacobi determinant, one has: 





a/^/2 






Jo 




'I 




Jo 




a" 






2 


J < 



ldYdX + 2 



a/V2 Jo 



h{X) 



1 dYdX 



^ (xj y A poo 

2 XdX + 2 h{X)dX 

Jo Jal\Pi 



h(X) dX, 



where 



X -Y X + Y 

u = = — , V = 



(5) 



V2 ' V2 ' 

and h{X) is determined by changing the variables u and v on the equation v = 
f{u) to the new variables X and Y, as follows. Since v = f{u) = sinh~^(coshu), 
then sinht; = coshu, so 



By making use of the addition and subtraction formulas for hyperbolic functions, 
as well as of the substitutions X' := and Y' :— Y/^/2, the above equality 

expands to 

sinhX' coshy' + sinhF' coshX' ^ coshX' coshF' - sinhX' sinhF' . 
By dividing both sides by coshX' coshF', one finds 

tanhX' + tanhr' = 1 - tanhX' • tanhF' , 
which can be written as 

1 - tanhX' 



tanhr' = 



1 + tanhX' 



(6) 



Now, for simplicity, let us define t tanhX'. After some algebra, Eq. ^ 
becomes 

1 -t- 

1 



Y' = -\n 
2 



l + t 



1 - 



1 - t 



= \nt. 
2 



which means that 



Y = h{X) 



V2 



In 



tanh 



i- 



By putting this function on the last integral of Eq. ([5|) , one has 



^2 a2 



In 



/\/2 



tanh 



X 
V2 



dX. 



The simple substitution z ~ X/^/2 promptly reduces this equation to 

/• oo 

/ In (tanhz) dz 

J a/2 



16 



(7) 

□ 



After this proof, we make the substitution t ~ tanh z in the integral in 
Eq. (O, above. When the result is expressed in terms of the dilogarithm func- 
tion, an interesting identity arises. 

Theorem 4 (Dilogarithm identity). Let us define Li2(2;) :— X^fcLi /k"^ , 
which is valid for all complex z with \z\ < 1 \S\. Then 

f r \ ( 1 \ In2(l + V2) 1 ^ 
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Proof. By applying the substitution t = tanh z [dt = sech^z dz) to the definite 
integral of Theorem [3l in the previous section, one finds 

"1 i„i rl 



a 
T 



TT 

16 



tanh(Q;/2) 



Int 



dt 



V2-1 1^ 



dt. 



since 



tanh 



tanh ^In \l \ + \/2 



a/2 



(8) 



\/2 - 1. 



By expanding the integrand in Eq. (|5]) in partial fractions, one has 



a 
T 



TT 

16 



1 



ln< 



1 r^-i Ini 



1 -i 



(9) 



These integrals can be reduced to special values of the dilogarithm function if 
one adopts the more general integral definition for this function namely 

In (1 - s) 



Li2(z) := - 



ds , 



/o s 

valid for z G C\[l,oo)|3 Note that, by expanding the logarithm in powers of 
z in this integral definition, one finds 1j\2{z) = /k"^! which agrees with 

the basic definition stated in the theorem for all complex z with \z\ < 1 0. By 
taking into account this summation form, it is clear that Li2(l) = 7r^/6, which 
implies that 



U2{z) 



In (1 - s) 



6 Ji s 
From this integral form, it is easy to deduce that 

Int 



ds . 



1 - t 



dt = Li2(l — z) 



and 



1 Int 

TTt 



dt 



-Li2(-z) -Inz In (-^ + 1) ^ ^ ■ 



By putting z — \/2 — 1 in these integrals and then substituting the resulting 
dilogarithms in Eq. one finds 



a 
T 



Tr- 
ie 



Li2 ( 1 - \/2 ) - Inf %/2 - 1 ) \nV2 



12 



Li2(2- \/2 



Li2(l-V2) _Li2(2-V2) _l , X 1,2-!^. (10) 

2 2 4 V y 24 ^ ' 



*We are considering here the principal branch of the dilogarithm function hi2{z), defined 
by taking the principal branch of In z, for which it has a cut along the negative real axis, with 
|arg(z)| < TT. This defines the principal branch of Li2(z) as a single- valued function in the 
complex plane cut along the real axis, from 1 to +00. 
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Now, by applying the Euler reflection formula for Li2(z) about z = i, which 
reads (see Eq. (1.11) in Ref. Q) 

Li2(l - z) = -Li2(z) -Inz ln(l - z) + y , 
to Li2(2 - %/2) in Eq. one finds 



Li9 



i2(V2-l)-Li2(l-%/2) =y-^ + ^-ln2(V2-l). (11) 
By applying the Landen's formula, namely (see Eq. (1.12) in Ref. Q) 

Li2(z) = -Li2(^^)-iln2 (l-z), 

which is valid for all complex z, z ^ [1, cx)), to Li2(l — V^), one has 

Liod- V2 ) = -Li,! 1- - In' 2. 



.12(1-^2) 



V2J 8 

Finally, by substituting this result on Eq. (ITT|) . one finds 

u.(v^-:).i,,(i--L)44-i„.. 



□ 



It should be mentioned that neither Mathematica (release 7) nor Maple 
(release 13) are able to generate any closed- form expression for this sum. Indeed, 
this identity is not found in standard bibliographic sources for dilogarithms [1,0, 
though it resembles the identity 

£(v2_i)+.(i--L).£, 

6 

where C{x) := — [Li2(a;) + i Ina; ln(l — a;)] is the normalized Roger's diloga- 

rithm, as found by Bytsko (1999) [13]. 

It seems plausible to extend the results obtained here in this paper to other 
known unit square integrals, such as J^J^ In {xy)/{l — x^y^) dx dy — —j ^(3) , 
J^J^lnx/{l-x^y^)dxdy = -|C(3), J^J^ Inixy)/ {1 - xy) dx dy - -2C(3), 

JqJ^ In (1 — xy) /(I — xy) dx dy = — C(3) , and many others (see, e.g., Ref. [HI). 
Similarly, it should also be interesting to include a third dimension in the hy- 
perbolic change of variables in order to investigate some triple integrals over the 
unit cube, such as 

m- dx dy dz — C(3) and III „ „ „ dx dy dz ~ - C(3) . 
l-xyz Jq Jf) Jq 1 - x^y^z^ 8 

Extensive work in both directions is already ongoing by the author. 
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Figures 




Figure 1: Mapping tlie unit square in the xy-plane (left-hand side) onto the hyperbohc 'quadri- 
lateral' in the uv-plane (right-hand side). The symmetry about the diagonal v = u, repre- 
sented by the inclined dashed line, is exploited in the text. Note that the hyperbolic region 
formed in the uv-plane is bounded by the coordinate axes and the curves v = f(u) and 
V = g{u). Both these functions, as well as the expressions for the change of variables indi- 
cated in the central part, are stated elsewhere in the text. According to Eq. ((Jjl, a evaluates 
to sinh-i(l) = \n(l + V2^ ^ 0.88137. 
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